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Viable modifications of gravity that may produce cosmic acceleration need to be screened in 
high-density regions such as the Solar System, where general relativity is well tested. Screening 
mechanisms also prevent strong anomalies in the large-scale structure and limit the constraints 
that can be inferred on these gravity models from cosmology. We find that by suppressing the 
contribution of the screened high-density regions in the matter power spectrum, allowing a greater 
contribution of unscreened low densities, modified gravity models can be more readily discriminated 
from the concordance cosmology. Moreover, by variation of density thresholds, degeneracies with 
other effects may be dealt with more adequately. Specializing to chameleon gravity as a worked 
example for screening in modified gravity, employing JV-body simulations of f(R) models and the 
halo model of chameleon theories, we demonstrate the effectiveness of this method. We find that a 
percent-level measurement of the clipped power at k < 0.3 h/ Mpc can yield constraints on chameleon 
models that are more stringent than what is inferred from Solar System tests or distance indicators 
in unscreened dwarf galaxies. Finally, we verify that our method is also applicable to the Vainshtein 
mechanism. 


Introduction Determining the nature of the accel¬ 
erated expansion of our Universe is a prime endeavor 
to cosmologists. In the conventional picture, the flat A 
cold dark matter (ACDM) concordance model based on 
general relativity (GR), a cosmological constant A con¬ 
tributes the bulk of the present energy density in the 
cosmos and drives the late-time acceleration. While al¬ 
ternatively, a modification of gravity may be responsible 
for cosmic acceleration, stringent limitations from experi¬ 
ments within our Solar System must be satisfied. A num¬ 
ber of screening mechanisms (J-Q have been identified 
that can suppress modifications of gravity in high-density 
regions to recover GR, while still generating significant 
modifications within lower densities on larger, cosmolog¬ 
ical scales. However, this suppression effect, along with 
other nonlinear effects, also prevents strong anomalies 
from manifesting in the averaged large-scale structure of 
our Universe ff| and limits the constraints that can be 
inferred on these gravity models from cosmology. 

Given the density dependence of the screening effect, 
in this Letter, we propose the downweighting of high- 
density regions in statistical observables such as the mat¬ 
ter power spectrum P(k) to enhance, or unscreen, the 
signatures of modified gravity and improve observational 
constraints. Such a weighting is conducted in the clip¬ 
ping method of Ref. [j], with the original motivation of 
facilitating the modeling of P(k) by reducing contribu¬ 
tions of high densities, where the assumptions of pertur¬ 
bation theory break down. As a worked example, we 
first focus on Hu-Sawicki @ f(R) gravity B, which em¬ 
ploys the chameleon screening mechanism |2|. We an¬ 
alyze effects on the power spectrum from clipping den¬ 
sity fields in numerical simulations of the model. Using 
the halo model of chameleon theories 0 , we then gen¬ 
eralize our findings to chameleon models with arbitrary 


gravitational coupling and exponents of the chameleon 
field potential. We also verify the applicability of clip¬ 
ping to Vainshtein screening jl|. 

Chameleon gravity .— We first specialize to the Hu- 
Sawicki f(R) (n = 1) model, where the nonlinear func¬ 
tion f{R) — — 2A — JroRq/R of the Ricci scalar R 
is added to the Einstein-Hilbert action. Here, bars 
denote quantities evaluated at the cosmological back¬ 
ground, zeros refer to present time, f R = df/dR is the 
additional, scalar degree of freedom of the model, and 
f R o = fn(z = 0). In the quasistatic approximation and 
for frto | <C 1, the modified Poisson equation becomes 
(see, e.g., Refs. 00) 

V 2 * = - ±5R{f R ), (1) 

where 6 denotes perturbations with respect to the cos¬ 
mological background, e.g., SR = R — R, and T = 
6goo/ (2goo)- The scalar field equation is given by 

V 2 Sf R = -^-Sp m + ^SR(f R ). ( 2 ) 

The chameleon mechanism works such that in a high- 
density region SR ~ 8nGSp m and hence Eq. © reduces 
to the Poisson equation of Newtonian gravity. In con¬ 
trast, at low densities SR ~ dR/df R \ R=R = 3m 2 Sf R , 
which when applied to Eqs. © and © in Fourier space 
yields the unscreened modified Poisson equation 


/c 2 'F = —47 tG 
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a 2 Sp m . (3) 


Hence, whereas in high-density regions gravity returns 
to Newtonian due to the chameleon mechanism, at low 
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densities and scales below m -1 , gravitational interactions 
remain enhanced by a factor of 4/3. In particular, Eq. ([3]) 
applies to the linear perturbation regime. Note that 
f(R) models correspond to a Brans-Dicke scalar-tensor 
theory with Brans-Dicke parameter to = 0, Jordan- 
frame scalar field ip = 1 + Jr, and scalar field potential 
U = (R/r — /)/2. More generally, the chameleon mecha¬ 
nism is realized for scalar-tensor models with scalar field 
potential U(ip)—A oc (1 — ip) a with a = n/ (n+1) £ (0,1), 
where the gravitational coupling is maximally enhanced 
by a factor of (4 + 2w)/(3 + 2w) with w > —3/2 jTcj. 
Importantly, although serving as very useful example for 
screening mechanisms, chameleon models do not yield a 
genuine self-acceleration of the cosmic expansion due to 

their gravitational modifications [13- 

In order to obtain accurate results in the nonlinear 
regime of the Hu-Sawicki f(R) model, we use dark mat¬ 
ter 7V-body simulations run in Ref. [l4j with the EC0SM0G 
code of Ref. EH, which uses particles and adaptive 
meshes to solve Eqs. m and ©• The background ex¬ 
pansion is taken to be equivalent to that of ACDM, ap¬ 
propriate for observationally interesting /rq values. We 
use simulations of the concordance model and f(R) grav¬ 
ity where \f R0 \ is 10“ 4 (F4), 10" 5 (F5), and 10” 6 (F6), 
all sharing an initial seed and cosmological parameters. 
Each simulation contains 512 3 particles in a box with 
L = 512/i _1 Mpc; h = 0.697, D m = 0.281, D b = 0.046, 
Da = 0.719, n s = 0.971, and amplitude of the matter 
power spectrum such that <ts = 0.82 in ACDM. Note that 
cj 8 is larger for f(R) gravity due to the enhanced forces 
and growth of structure. An initial power spectrum for 
the simulations was generated using MPGRAFIC [l6|. The 
particle mass in each case is ~ 7.80 x 10 10 7?.“ 1 M 0 . Each 
simulation has exactly the same initial power spectrum 
( z\ = 49) and differences between models are confined 
to different strengths of enhanced perturbation growth 
at late times and different strengths of screening. For 
the extrapolation of nonlinear physics from f(R) grav¬ 
ity to more general chameleon models with u / 0, we 
employ the halo model of chameleon theories developed 
in Ref. [Ic|, which accounts for the density dependence 
of the effective gravitational coupling in the nonlinear 
regime and provides matter power spectra that are in 
good agreement with measurements in f(R) IV-body sim¬ 
ulations (see Fig. 4 in Ref. E3). 

Clipping the density fields. — The spatial distribution 
of matter on cosmological scales may be quantified by 
the fractional overdensity field S(x) = pm/fim — 1- We 
construct the density fields from the simulations using 
a cloud-in-cell interpolation on a 256 3 Cartesian mesh in 
each cell. Clipping is a local density transformation char¬ 
acterized by enforcing a maximum fractional overdensity 
So such that [7] 


S c (x) 


S 0 , S(x) > S 0 , 
S(x), S(x) < So- 


(4) 


We shall also make use of applying a minimum instead 
of a maximum threshold, which is equivalent to clipping 
the negative field —S(x). While this may prove more 
challenging to apply to real data, due to the lower signal 
to noise associated with cosmic voids, it will serve as a 
useful validation test in our simulations and support for 
the concept of weighting to unscreen or screen gravita¬ 
tional modifications. We quantify the clipping strength 
in terms of the fractional loss of power in the lowest k- 
bin, applying a simple iterative procedure to determine 
the threshold <5 0 required to establish the desired frac¬ 
tion. Defining clipping strength directly in terms of So 
is another possibility, but complicates the comparison of 
results from fields with a different choice of smoothing 
length. 

Given the nature of the chameleon mechanism, the 
clipping transformation promises the extraction of more 
information on the gravitational physics, as it allows us 
to focus on the less dense, unscreened regions of the Uni¬ 
verse, where there exists a greater difference from GR. In 
the left panel of Fig. [I] we compare the deviation in power 
between the F6 and ACDM simulations as a function of 
clipping strength. With increased clipping strength, the 
difference in P{k) is enhanced. At k = 0.5 h /Mpc this 
difference more than doubles, greatly facilitating the dis¬ 
crimination between the two models. In practice, the 
amount of clipping that can realistically be applied will 
depend on the noise levels, limited by the number density 
of galaxies. To see the reverse of this effect, we also apply 
the clipping transformation to the negative field — <5(x). 
Since this is clipping voids, this reduces the contributions 
from unscreened regions to P{k). Therefore, as expected, 
stronger clipping leads to convergence between the mod¬ 
ified gravity and ACDM density fields. Note that we do 
not expect a full unscreening of modified dynamics by 
clipping high-density peaks since this only affects self¬ 
screening but not environmental screening 53- In the 
opposite case of clipping low-density troughs, however, 
only self-screened regions with GR dynamics remain. 

In the middle panel of Fig. [TJ we show the ratios of 
the power spectra of F4-6 to that of ACDM simulations 
and the corresponding ratios after clipping has been ap¬ 
plied such that the large-scale power is reduced by 50%. 
We also show the results of instead applying a logarith¬ 
mic density transform 17]. Both local transformations 
suppress contributions from the densest regions of the 
field. Note that while the shape recovered from the log¬ 
arithmic transform is closer to linear theory, for sparsely 
sampled fields the logarithmic transform becomes unsta¬ 
ble whereas cli ppin g is insensitive to the number den¬ 
sity of sourc es [1811 . Regarding the proximity to linear 
theory, Ref. |18[ demonstrated that the clipped ACDM 
matter power spectrum is well described by a linear com¬ 
bination of the linear and one-loop contributions, where 
higher-order terms are strongly suppressed. Thereby, in¬ 
creased clipping downweights the one-loop relative to the 
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FIG. 1: Left: Fractional difference between simulated matter power spectra of f(R) gravity (|/flo| = 1CP 6 ) and ACDM at 
two different fc-bins. By clipping high densities, contributions from self-screened regions are downweighted and unscreened 
features enhanced (dashed curves). In the reverse case of applying a minimum threshold to the density field, contributions from 
unscreened low densities are downweighted, enhancing screened regions, and the power spectra converge (dotted lines). Middle: 
Fractional departure from the ACDM power spectrum in three different f(R) models for unclipped fields (solid lines), linear 
perturbation theory (dashed lines), clipped fields with power at the largest scales reduced by 50% (dot-dashed curves), and 
logarithmic density transformation with matching to the ACDM power at the lowest fc-bin (dotted curves). Right: Clipped f(R) 
power spectra from simulations (circles) compared to weighting linear and one-loop contributions (dashed curves). Increased 
clipping strength enhances the weight of the unscreened linear (dotted curves) relative to the one-loop power. 


linear contribution. We verify the recovery of linear the¬ 
ory for f(R) gravity in the right-hand panel of Fig. |I] 
using F4 simulations, where chameleon screening affects 
P[k) at k > 0.1 h/Mpc (see Fig. 2 in Ref. [(|). With 
increased clipping, P(k) more closely reflects linear, un¬ 
screened, theory with a reduced relative weight of the 
one-loop correction. 

Differential clipping to break degeneracies .— Effects on 
the total matter power spectrum from f(R) modifications 
of gravity, and hence other chameleon models, have been 
shown to yield some degeneracy with effects from varying 
the total neutrino mass [l{J [io( or from baryonic feed¬ 
back processes [2l|. It is important to note that while 
degenerate in P{k ), these contributions are not affecting 
a particular range of densities in the same manner as the 
chameleon modification. More precisely, nonvanishing 
neutrino masses suppress P(k ) predominantly linearly, 
i.e., with negligible preference on density, and baryonic 
feedback processes emanate from high-density regions. 
Hence, a combined analysis of clipped power spectra em¬ 
ploying different density thresholds can break degenera¬ 
cies between the different contributions. Furthermore, in 
a parameter estimation analysis of chameleon models, we 
need to allow for a variation of cosmological parameters. 
In particular the growth of structure can be strongly al¬ 
tered by variations of a 8 or H m . Using HALOFIT [22^ . 
we estimate effects of varying these parameters on the 
clipped power by comparing signatures in the linear and 
nonlinear P(k). While enhancing the amplitude of P(k) 
with increasing as, the shape of the nonlinear enhance¬ 
ment resembles that of a chameleon model. Hence, with 
absent information on the absolute amplitude of P(k ) 
due to galaxy bias, there clearly is a degeneracy between 


variations in as and /rq. However, the ACDM case with 
a larger value of as, will experience a strong reduction 
of the enhancement in power at large k after clipping. 
In contrast, the f(R) power spectrum increases the en¬ 
hancement at large k after clipping, becoming more lin¬ 
ear and unscreened. Therefore variations in /rq and as 
respond to clipping in a qualitatively different manner, 
breaking the degeneracy. Modifications in the shape of 
P(k) due to changes in fl m are qualitatively different 
from f(R) modifications, e.g., changing baryon acous¬ 
tic oscillation features. Although a partial suppression 
in the change of power attributed to variations is 
seen in nonlinear compared to linear theory, it does not 
reproduce the strong screening effect of chameleon mod¬ 
els. Finally, note that since focusing on differences in 
the shape of P{k), effects of linear galaxy bias can be 
neglected. However, through redshift-space distortions, 
in chameleon models the ratio between galaxy and dark 
matter density becomes scale dependent. Since adding to 
the deviations between the shape of modified and ACDM 
galaxy power spectra [23], we conservatively assume this 
ratio to be constant when estimating potential observa¬ 
tional bounds on chameleon models. 

Outperforming Solar System constraints .— The re¬ 
quirement that the Milky Way dark matter halo screens 
the Solar System sets a constraint on the chameleon field 
amplitude of \(po — 1| ^5x 10 -6 / (6 + 4w) [Ic|. Similarly 
strong constraints can be obtained from the absence of 
deviations in luminosity distances from different types 
of distance indicators in unscreened dwarf galaxies [24j . 
Current cosmological constraints are about 2 orders of 
magnitude weaker than Solar System bounds [25l - !27| . We 
refer to Ref. (28| for a review of constraints on chameleon 
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gravity. Having analyzed the effect of clipping on the 
matter power spectrum, we estimate the constraints that 
can be inferred from applying this method to observa¬ 
tions of galaxy clustering. While fractional errors in the 
measurement can be kept approximately constant, clip¬ 
ping reduces systematic errors from modeling uncertain¬ 
ties of the nonlinear structure contributing at large k 
such that constraints on cosmological parameters can be 
improved (see Ref. ;2f|). We assume that a future mea¬ 
surement of the clipped P(k) can discriminate 1% devi¬ 
ations from the fiducial ACDM shape at k < 0.3 h/ Mpc. 
Note that this is a conservative estimate with surveys 
already yielding subpercent-level measurements of the 
acoustic features §3(|. We compute the modified nonlin¬ 
ear power employing the halo model of chameleon the¬ 
ories [lOj. As clipping chameleon densities mainly re¬ 
moves regions where GR is recovered, we approximate 
the clipped chameleon power spectrum by removing the 
difference of undipped to clipped power obtained from 
the ACDM simulation. For k < 0.3 /i/Mpc, this simple 
approach reproduces the absolute clipped power for F4 
and F6 at the few percent and permille level, respectively, 
and recovers the measured fractional difference to ACDM 
within 20%, increasingly underestimating it from F4 to 
F6. Employing this method and varying the background 
field value ipo, the coupling strength set by w, and the ex¬ 
ponent of the scalar field potential a, we set constraints 
where the deviation between the clipped chameleon and 
ACDM P(k ) at k = 0.3 h /Mpc exceeds 1%. 

We present our results in Fig. [2j comparing them to 
Solar System, astrophysical distance indicator, and cur¬ 
rent cosmological bounds. We find that for the simulated 
Hu-Sawicki model (a = 0.5), clipping constraints on /rq 
can improve upon existing cosmological constraints, us¬ 
ing clusters 251, the matter power spectrum j2(|, or cos¬ 
mic microwave background with redsliift-space distortion 
data [ItJ , by 2-3 orders of magnitude and outperform the 
Solar System and astrophysical bounds as well as local 
tests of the equivalence principle [2j, |H, 31). For smaller 


values of w, corresponding to larger force modifications, 
the improvement of constraints from clipping over Solar 
System bounds is even greater. This is not surprising 
since larger force modifications also imply a more effi¬ 
cient screening of the Solar System region whereas gravi¬ 
tational dynamics in unscreened low densities is modified 
even more strongly. Constraints in this region of parame¬ 
ter space may also confirm or rule out chameleon models 
as an explanation of the observed cored density profiles 
of dwarf spheroidal galaxies in the Milky Way [32|. In 
the opposite limit of increasing to, i.e., weakening grav¬ 
itational coupling, Solar System bounds strengthen and 
surpass the decreasing clipping constraints. Importantly, 
these power spectrum constraints clearly depend on the 
exponent of the scalar field potential a. 

Vainshtein screening A similar density dependence 
to the chameleon mechanism enters the effective gravita- 



FIG. 2: Forecast of upper bounds on the chameleon field am¬ 
plitude 1 — tpo from the clipped matter power spectrum as a 
function of gravitational coupling strength set by the Brans- 
Dicke parameter to and for different values of the exponent of 
the chameleon field potential a. The constraints are compared 
to upper bounds inferred from Solar System [32}, astrophys¬ 
ical distance indicator 24}, and cosmological (cluster) obser¬ 
vations [25} (cf. Ref. [28|]). For a range of chameleon model 
parameters, clipping can provide the strongest constraints. 


tional coupling at nonlinear scales in models exhibiting 
a Vainshtein jjj screening effect §33}. While the modifi¬ 
cations of P(k) are comparable, the Vainshtein mecha¬ 
nism originates from derivative self-interactions in con¬ 
trast to the scalar field potential in chameleon gravity, 
causing some differences in the screening behavior [341 ]. 
Chameleon screening depends on the gap between an ob¬ 
ject’s internal and exterior scalar field values set by the 
corresponding densities. Hence, an object can be self- 
screened or environmentally screened by a large local or 
environmental density, respectively. In contrast, in the 
Vainshtein mechanism, screening depends on morphology 
of the source but not on its environment, introducing a 
large screening scale around an object, suppressing modi¬ 
fications to its own gravitational field but still responding 
to external field modifications. We verify that the clip¬ 
ping method is also applicable to the Vainshtein mecha¬ 
nism by employing 50% clipping to the nDGP §35} sim¬ 
ulations of Ref. [34J . We find an enhancement of A P/P 
(cf. Fig. []} of 70%, 60%, and 20% at screened scales of 
k ~ 1.5 h /Mpc for models with cr 8 values equal to the F4, 
F5, and F6 f(R) scenarios, respectively. We also expect 
our method to be applicable to further screening mecha¬ 
nisms like the symmetron §4] or k-mouflage [3], however, 
not to linear shielding [§]. 

Conclusions .— Modifications of gravity potentially ex¬ 
plaining cosmic acceleration need to employ a screening 
mechanism that allows modifications in low densities at 
large scales while suppressing them in high-density re¬ 
gions such as the Solar System. Such screening mech¬ 
anisms also suppress anomalies in the averaged cosmo¬ 
logical large-scale structure, limiting observational con¬ 
straints inferred on such models. By clipping the density 
fields at a maximal threshold, contributions of screened 
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high-density regions to the matter power spectrum can 
be downweighted, enhancing modified gravity effects. We 
applied this method to chameleon models with particu¬ 
lar emphasis on f(R) gravity to demonstrate that it can 
improve cosmological constraints on the models to a level 
stronger than the currently most stringent bounds from 
Solar System and astrophysical tests. We also verified 
that clipping is applicable to Vainshtein screening and 
expect it to be employable for further nonlinear screen¬ 
ing mechanisms. 
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